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Abstract 

Hamiltonian tridiagonal matrices characterized by multi-fractal spectral mea- 
sures in the family of Iterated Function Systems can be constructed by a recursive 
technique here described. We prove that these Hamiltonians are almost-periodic. 
They are suited to describe quantum lattice systems with nearest neighbours cou- 
pling, as well as chains of linear classical oscillators, and electrical transmission 
lines. 

We investigate numerically and theoretically the time dynamics of the systems 
so constructed. We derive a relation linking the long-time, power-law behaviour 
of the moments of the position operator, expressed by a scaling function /3 of the 
moment order a, and spectral multi-fractal dimensions, Dg, via /3(a) = Di-a- We 
show cases in which this relation is exact, and cases where it is only approximate, 
unveiling the reasons for the discrepancies. 
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1 Introduction 

Usually, the study of almost/quasi-periodic systems starts by assigning a suitable rule 
for building a quantum Hamiltonian operator, and then proceeds to the determination 
of its spectral quantities (which is frequently a hard task) and of the time dynamics it 
generates. In so doing, multi-fractal energy spectra have been frequently observed, and 
anomalous characteristics of the time evolution have been exhibited |, ||, |], 

These findings raise the question if multi-fractal spectra are typical in almost/quasi- 
periodic systems and, vice versa, if almost/quasi-periodicity is always associated 
with singular continuous spectral measures. The relations between this pair (spectral 
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multi-fractality and Hamiltonian almost-periodicity) and the time dynamics generated 
via Schrodinger's equation are also interesting, and intricate: do the former always imply 
anomalous scaling relations of physical quantities like, for instance, the expectation value 
of the position operator ? Can we make any quantitative statement to this effect ? 

In this paper, we employ a new algorithm for deriving a Hamiltonian operator (called 
a Jacobi matrix because of its mathematical nature) with a pre-assigned spectral measure 
in the vast class of Iterated Function Systems (I.F.S.). This technique provides us with an 
ideal patient for our surgical table, who can be fully dissected and analyzed. In particular, 
we provide evidence that Jacobi matrices associated with I.F.S. are almost periodic, and 
we argue that this is likely to be the typical case in a large class of measures with fractal 
support. 

The analysis which can be carried out in this example permits us to compute exactly 
the asymptotic behaviour of the wave-function projections, for short and long times. By 
introducing a renormalization approach in the theory of orthogonal polynomials, we also 
derive a relation linking the asymptotic power-law growth of the moments of the position 
operator and multi-fractal generalized dimensions. This theory explains the phenomenon 
that we have termed quantum intermittency. 

The specific properties of I.F.S. are crucial for our theory; Yet, since particular I.F.S. 's 
can be found so to approximate arbitrarily well (in a technical sense) any "fractal" 
measure HSl. 



J, ^, |10|, |Tl|], the results obtained in the I.F.S. class may have a much wider 
generality. 

We shall present our results as follows: in Section II we introduce the general formal- 
ism of I.F.S. and of Jacobi matrices, employed to solve the inverse problem of finding 
an Hamiltonian with a given spectrum. This formalism is then applied in Sect. HI to 
derive a stable solution algorithm. The almost-periodic properties of the Hamiltonian so 
determined are studied numerically in Sect. IV, and the intermittent quantum dynamics 
it generates is then discussed in Sect. V and VI. The Conclusions summarize the work 
and present some previews on further investigations. 



2 I.F.S. and their Jacobi Matrices 



Systems of linear iterated functions [|T2], 0, [T^, [T^ are finite collections of maps 

<j)i{x) := 6iX + /3i, i = l,...,M, (1) 

where 5i, /3j are real constants, and where the contraction rates Si have modulus less than 
one. Without loss of generality, we may assume that each 0j maps [0, 1] into itself, and 
that 01 (0) = 0. 

A probability, tTj, is associated with each map: tTj > 0, = 1- Employing these 

probabilities, a measure over [0, 1] can be defined as the unique positive measure satisfying 
the balance property 



^0 Jo 



(2) 



for any continuous function /. This measure is supported on A, the subset of [0, 1] which 
solves the equation 

A= U <I'^(A)^ (3) 



i=l,...,M 
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The set A is invariant under the action of shrinking it to smaller copies of itself, and 
glueing them together. Because of eq. the geometry of this set is typically fractal 
(except for special choices of the map parameters); In turn, the balance relation (|^) is 
responsible for the multi-fractal properties of the measure /i. In fact, let us consider a 
disconnected I.F.S., that is to say, one for which the sets (pi{A) do not intersect each other. 
Under these circumstances, the multi-fractal properties of the balanced measure are easily 
computable: the spectrum of generalized dimensions Dq follows from the equation 



M 



E-ir = i' (4) 



whose unique real solution defines r as a function of q, and leads to Dq = ^3^. In virtue 
of this relation, one can tune the map parameters to obtain various multi-fractal spectra. 

The problem of determining a Hamiltonian possessing /i as spectral measure can be 
solved [|16] considering the set of associated orthonormal polynomials, {pn}'- 



Pi{x)pk{x) d^{x) = 6i^k- (5) 

In fact, any such set of polynomials is characterized by a three-terms recurrence relation 
which can be written 



xpj (x) = Tj+ipj+i (x) + AjPj (x) + rjPj_i{x) , (6) 

or, in matrix form 

Hp{x) = xp{x). (7) 

In the above, p{x) is the vector whose components are the orthonormal polynomials 
evaluated at site x, and H is the Jacobi matrix, which is constructed as the real, sym- 
metric, tridiagonal matrix whose diagonal and outer diagonals are the vectors Aj and rj, 
respectively: 

Hi^i = Ai, -ffj+ii = -ffii+i = Tj+i, z = 0, 1, (8) 

H defines a nearest neighbours lattice system, with site energies Ai and hopping constants 
Tj. Similarly, H can describe a linear array of masses coupled by springs, and also an 
electrical transmission line, whose characteristics vary from one element to the next. 

Standard theory proves that, letting the Jacobi matrix H act in I2 (the space of square 
summable sequences, whose canonical basis will be indicated by {cq, ei, . . .}) the spectral 
measure of H with respect to the vector cq (the local density of states of physical jargon) 
is precisely fj,: in fact, one has 

{eo,g{H)eo) = J g{x) dfi{x), (9) 

for well-behaved functions g. This is the theoretical solution of the inverse problem we 
have proposed. In order to translate it into a practical solution, we need to compute the 
Jacobi matrix coefficients starting from the measure i.e. from the map parameters 
defining the I.F.S. 
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3 A Stable Technique for Computing I.F.S. Jacobi 
Matrices 



The problem of constructing the Jacobi matrix associated with I.F.S. measures is hard, 



and the usual techniques oi polynomial sampling |17, IS, IS] are plagued by exponentially 



increasing errors which allow only computation of very few Jacobi matrix coefficients 



20]. Alternatively, the sole technique available so far has been an algebraic procedure 
programmed in MAPLE hj Vrscay [^. Yet, it is severely limited by memory and time 
requirements. To overcome these difficulties we have devised a direct algorithm applicable 
to I.F.S. measures. 

We first observe that, for any n, 

n 

PniMx))=T.^llPlix), 2 = 1,...,M. (10) 
1=0 

This is immediate, since pn{(f)i{x)) is an n-th degree polynomials which can be expanded 
on the first n orthogonal polynomials. Less immediate is to derive a recursive rule for the 
coefficients F";, I = 0, . . . ,n. It turns out that, at fixed n, they can be determined from 
the map parameters, and from the Jacobi matrix entries Aj, for j = 0,1, . . . ,n — 1, and 
Tm, for m = 0,1, . . . ,n. In fact, (dropping for simplicity the map index i) we have that 
Po{6x + P) = po{x), and hence Fq = 1. Suppose now that Ff is known for = 0, . . . , n — 1 
and all relative /'s: from eq. (^ we obtain the complete decomposition of Pn{^x + /9) 
over pi, I = 1, . . . ,n\ 

rnPniSx + P)= {13- A„_i) EITo' ^?-'pi{x)+ 

+5J:?=o rr\ri+iPi+i{x) + Aipiix) + ripi^i{x))+ (11) 

-r^-iE?=o'rr'piix). 

Equation (|TT|) allows now the determination of the coefficients F". 

We observe that the highest order polynomial, Pn, appears twice in the above equation, 
always in the form of the product r„p„: hence, the coefficients in the expansion of the 
polynomial r„p„ can be determined without knowing r„. Therefore, if we let Pn{x) = 
fnPn{x), a second decomposition can be written as 

n-l _ 

P„(0.(X)) = pn{x) + ni P^(^)^ (12) 

1=0 

where the coefficients F can be computed recursively from eq. ([TT|), on the basis of the 



knowledge of only Aj, rj, for j = 0, 1, . . . , n — 1. 

We can now compute the non-diagonal entries of the Jacobi matrix: from eq. (^ we 
write 

^1= Pn{x)xpn^i{x) dfl. (13) 



Hence, using the balance property (^ and eqs. ([10|) , (p!2|) this becomes 

AI „ n—1 n—1 n—1 



^2 



(s.x + A)[E E ri^rir'pm{x)pi{x) + f ^„F^r'pn(x)pK^)]rf/i- (i4) 



i=l m=0 1=0 1=0 
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Again, we can use the recurrence relations (^), to get 

AI 

Y^-Ki (Bi + Q + D,), (15) 



n 

i=l 



where we have put: 



n-1 



= + S^Al)Tl,TlT\ (16) 

n-2 

a = 5. E n+i(f^,r]:^ + f^,+iri;r'), (ir) 



/=0 



and 



-^j — y«J- j,n-^ i,n-l' n- 



Because of contractivity of the maps, < 1- Therefore, (and hence r„ > 0) 

can be computed from eq. (0), on the basis of the knowledge of the coefficients in the 
expansions (|l^) of order n — 1, of order n in ([l^), of the map parameters, and of the 
matrix entries Aj, rj, for j = 0, 1, . . . , n — 1. 

A similar trick allows the computation of the diagonal entries An] We use eqs. @ 
and (BI) (integrals are taken with respect to fi): 



. M . M . n 

An = J Xpl{x) =J2^i J {SiX+(3i)pl{diX+f3i) = ^T^i J (SiX+Pi) J2 ^i,l^lmPlix)Pm{x) ■ 
i=l i=l m,/=0 



(19) 

Using the orthonormality properties of the sequence Pn, and the recurrence relation, eq. 
(D, we get 



M n n-1 



An = E^aE (r-n.)' (A + ^.^n.) + E nm^m^Mrn. + r„.+l)], (20) 



1=1 m=0 m=0 



thereby determining An as a function of the coefficients in eq. (p!0[) of order n fixed, of 
the map parameters, and of the matrix entries Aj, for j = 0,l,...,n — 1, and r^, for 
m = 0, 1, . . . , n. 

These results can be properly chained into an iterative construction of the Jacobi 
matrix H: The algorithm is structured as follows: 

• Initialization. At the first step, we have Aq = fii, tq = 0, Fq = 1. The first order 
moment of /i, /ii, can be simply computed from eq. (H). 

• Iteration. Suppose that Ai, ri, and F' are known for / = 0, 1, . . . , n — 1. Then we: 



Compute F". We use equations (IT^ - p^ . 
Compute Tn- We use eqs. (|I^ - |TS|). 
Compute F". This is immediate at this stage. 

Compute An. We use eqs. (|T^ - ^). Then we iterate the procedure. 



5 



Graphically: 



TTl— 1 




fO, ■ ■ ■ 

Ao,...,A, 



:2i] 



In a separate work |^ we have analyzed the reasons of the failure of classical poly- 
nomial sampling ||T7| , |19[ when applied to singular measures, and assessed the numerical 
stability of the recursive algorithm presented above. We have observed a polynomial 
error propagation with respect to matrix order for the recursive algorithm, while using 
the classical algorithms the error growth was found to be exponential. 



4 Almost Periodicity of I.F.S. Jacobi Matrices 

Having devised a stable solution of the Hamiltonian inverse problem, we can study the 
properties of large Jacobi matrices. Fig. 1 shows an I.F.S. measure, one of its orthogonal 
polynomials, and the begining of the sequence of r„ coefficients. Let us focus our attention 
on the last. 

We can clearly observe a zero frequency (the average value) , a tt frequency (flipping up 
and down), and clearly other frequencies are present in the sequence. A Fourier analysis 
is simply effected writing 

rn = Y.Fke'"^'- (22) 

k 

This sum may not converge in the usual sense, and it might have to be replaced by an 
integral in the case of a continuous component in the "spectrum" of the sequence r„. If 
the continuous component is absent, the system is almost-periodic. Within this case, if 
the set of frequencies u>k can be derived from a flnite set of periods, the sequence r„ is 
quasi-periodic: that is, this is the case if there exist suitable Qi, . . . ,Qp such that for all 
k the frequency cj^ can be written cu^ = niQi + . . . + Upflp, for integer rii, . . . ,np. 

A numerical, fast Fourier analysis of the sequence r„ is presented in Fig. 5, where 
peaks in the distribution of with a clear hierarchical structure are observed. These 
peaks seem to suggest the presence of a point component in the spectrum of this sequence. 
Yet, care has always to be exerted to assess this fact numerically. To obtain a further 
piece of evidence we performed an analysis of the phase of F^. around these peaks, like 
that shown in Fig. 2, and found a vr discontinuity, which indicates |^ that they are 



indeed related to a point component. The sequence r„ is therefore almost periodic. 

Since no simple rational relation among the peak sequences seems to hold, numerical 
evidence seems to suggest that the sequence is not quasi-periodic. Our numerical inves- 
tigations have shown that these characteristics are typical in the class of Hamiltonian 
associated with I.F.S. measures, supported on Cantor sets. In view of the approximation 
properties of I.F.S. measures, this result is likely to be much more general: indeed, in 



the family of Jacobi matrices associated with real Julia sets |23[,[^], which can be well 
approximated by I.F.S., limit periodicity of the sequence r„ has been proven directly |25 
The problem of a formal proof is therefore open. 
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5 Quantum Dynamics of Almost Periodic Lattice Sys- 
tems 



Jacobi matrices generate a quantum dynamics in I2 via Schrodinger's equation, 

t^ = H^^, ^(0) = eo:=(l,0,...). (23) 

The initial state of the evolution, cq, is the zeroth lattice state. In oscillator terms, this 
corresponds to a situation where the first mass is displaced from its equilibrium position, 
while all the other masses are at rest in their equilibria. In electrical terms, the current 
(or the voltage) is non-zero only in the first element of the transmission line described 
by the Jacobi matrix H. 



The solution of Schrodinger equation can be formally obtained as p6|, ^ 

Cn{t) := (e„, e-**^eo) = / e-**>„(x) dfi{x), (24) 



where c„(t) is the component of ■0(t) at the n-th lattice state. Equation (|2^ ) shows that 
this component is the Fourier transform of the orthogonal polynomial pn with respect to 
the spectral measure /z. This fact allows us to derive important results. 

Firstly, the asymptotic behaviour for small t can be controlled as follows: |c„(t)P ~ 
^2". In fact, 

I 00 / 



c„(t) = J dfi{x) p„(x) ^ ^jr~^^ = £ / dn{x)p„{x)xK (25) 

1=0 ' l=n 



Because of the orthogonality properties of the set p„ this expansion begins with / = 
which proves the result. 

Secondly, in the infinite time limit, denoting by SniT) the time average of |cnp up to 
time T, 

^n{T) = ^j\n\\t) dt, 

we have that 

SniT) ~ T-^^ (26) 

for all n, a result which involves the correlation dimension D2 of the fractal measure 
Ij,. The case with n = is implicitly contained in Bessis et al. and was originally 
proposed in the present context by Ketzmerick et al. 0. Successively, it has attracted a 
lot of attention, mainly from cultors of mathematical rigour. Our generalization has the 



advantage of requiring a simple proof, via the usage of the Mellin transform, as in |28 
In fact, we write 

SniT) = J d^ix) J d^iy) Pnix)Pniy). (27) 

To find the asymptotic behaviour of eq. (p7|), we take the Mellin transform, M„(2;), of 
^niT): 



Mniz) = J T'-'SniT)dT = Gx J dflix) J dM ^]^f_^y\^!^ = Giz) x Eniz), (28) 
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where G{z) = T(z — l)sin(|(z — 1)), and where En{z) is defined imphcitly by the last 
equahty. The dominating power law in the long time behaviour of Sn is determined by 
the divergence abscissa of Mn{z): that is to say, S'„(T) ~ T~^, where w is the largest 
real z for which Mn{z), hence En{z) converges. It is apparent from eq. ( pH]) that the 
divergence of En is piloted by the small scale structure of the measure /i. Because the 
polynomials Pn are smooth functions, with bounded derivatives on the support of /x, the 
divergence abscissa of En is the same for all n, and, in particular, it coincides with that 
of £^0- Eq[z) is known as the generalized electrostatic energy of the measure n and its 
divergence abscissa is known to be D2 [^, the correlation dimension of the measure //. 

It is important to remark that the domains of validity of the asymptotic expansions 
just derived are not uniform in n. This adds to the difficulty of the problem to be 
discussed in the next Section. 



6 Renormalization Theory of Quantum Intermittency 

An important characteristics of the quantum motion introduced in the previous section 
is the way it spreads over the I2 lattice basis, {e„}. In fact, in oscillator terms, spreading 
corresponds to energy transmission along the linear chain, be it mechanical or electrical. 
In quantum mechanical terms, it corresponds to unbounded motion of the lattice particle, 
of the kind treated only qualitatively by R.A.G.E. theorems. To gauge this phenomenon, 
we define the moments of the position operator h: 

v^{t) ■= m),h^i,{t))=Y,n'^\cn{t)\\ (29) 

n 

Their asymptotic behaviour follows a power law, 

Vait) ~ t"^ (30) 

where /? is a non-trivial function of the moment order a. In 27| we found that 13 



is convex, non-decreasing, and non-constant even in the case of a one-scale Cantor set, 
characterized by trivial thermodynamics: this is what we call quantum intermittency. 
Corrections to eq. ( PDD can also be observed in the form of log-periodic oscillations of 
fQ,(t), super-imposed to its leading behavior. They can be explained by the Melhn-type 
analysis presented in the previous section. 

We can estimate the function /3(a) on the basis of simple renormalization group 
considerations. For simplicity, let us consider an I.F.S. with M maps, of equal probability 
TTj = -j^. Let this I.F.S. be non-overlapping. Then, let / be the smallest interval containing 
A, the I.F.S. attractor, and let J; be the image of / under the map 0;. Clearly, J; fl/m = 
if / 7^ m, and the measure restricted to /; is a linearly rescaled copy of the original. 
Then, as a first approximation, we can assume that the orthogonal polynomials of the 
restricted measure are also obtained by linear rescaling of the original polynomials: 

Mn 

VMnm^)) = E Kk^Pkix) ^ <P„(X), (31) 
fc=0 



where erf = ±1. In other words, we assume a very simple form for the coefficients F 
which amounts to making a renormalization ansatz. 



Mn 
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Let us now consider SMn{T), as defined above. Because of tlie balance property (0), 
it can be written 

SMn{T) = '^l^rn I d^{x) f dn{y) ^ ■ ^^7^ ^ P^n {4>l (x) )pMn{(t>m {v) ) • 

it^Li J J T{cf)i{x) - cf)m{y)) 

(32) 

In tfie previous equation, 4>i{x) and 4>m{y) belong to /; and /„, respectively. If / 7^ m, 
these intervals are separated by a finite gap. As T tends to infinity, these contributions 
tend to zero as T~^. We can therefore retain only the diagonal terms in eq. (0). 

If we now employ the approximate estimate (|3lD in the r.h.s. of eq. (p2D we can write 

M 

SMn{T)=J27rfSr,{6iT). (33) 
1=1 

This too is a sort of renormalization equation which links the wave-function component 
at site Mn and time T to the component at site n and at shorter times 5iT. When 
inserted in eqs. (|29|) , eq. (|33|) implies that the growth exponent P associated with the 



averaged moments Va via eq. (0) must satisfy the relation 

M 

1 = M"-i Y.^?^- (34) 
1=1 

Comparing this result with eq. (^) we obtain the crucial equation 

f3{a) = (35) 



which links multi-fractal properties and time dynamics. In particular, eq. ( p5|) implies 
that P{0) = Di, which is consistent with the rigorous result j3{0) > Di 0. Notice 
that /3(0) can be defined by a limiting procedure on /?(«), or by the evolution of the 
logarithmic moment. We have also (3{1) = Dq. 

Because of the rough approximation involved in eq. (|3T|) , and because for the validity 
of eq. ( P^ ) both cmu and c„ need to be in their asymptotic regimes, we do not expect eq. 
(p5D to be always exact. Indeed, in Fig. 4 we have considered a family of I.F.S. measures. 



characterized by M = 2, (^2 = |, /?i = 0, = |, tti = |, and = |. The contraction 
rate 61 is allowed to vary in the range [.2, .5], which implies a significant variation both in 
the structure of the support of the balanced measure and in its multi-fractal properties. 
Plotted in Fig. 4 are the scaling exponents (3{0) and /?(!), compared with the multi- 
fractal dimensions Di and Dq, respectively. We observe a substantial agreement between 
the two data sets, dynamical and multi-fractal. Numerically, the discrepancy is always 
less than five percent. We can therefore conclude that the relation ( ^51) catches some 
essential part of the physics. Yet, the situation is more complicated, as the following pair 
of examples show. 

The first is a magnificent counter-example. Let us consider a new class of I.F.S. 
measures (and related Hamiltonians) characterized by M = 2 and by a particular choice 
of the weights: 

TT, =5f, J = 1,2 (36) 

where D is the (constant) value i^^^§z^^- This choice originates what is called a uniform 
Gibbs measure. The first of such I.F.S. is that of Figs. 1 to 3, and D is its fractal 
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dimension. Indeed, all I.F.S. with the property (^) are characterized by the same flat 
thermodynamic function Dg = D. Clearly, because of eq. (0), and because tti + 7r2 = 
1, only one parameter among the map weights and contraction rates is left free. By 
varying this parameter we can construct different I.F.S. measures, with the same flat 
thermodynamics. What are then the corresponding dynamical exponents I3{a) ? The 
approximate relation p5| ) predicts /5(a) ~ D for all a. 

In Fig. 5 we have considered: a: The I.F.S. with 5i = 82 = |, tti = 712 = |, which is a 
"pure" Cantor Set. b: The I.F.S. with 5i = .5090, 62 = .2978, and tti = |. c: The I.F.S. 
with 61 = .5293, 62 = .2802, and tti = .6180. d: The I.F.S. with Si = .6033, 62 = .2196, 
and TTi = .6823. The first observation we can draw from this figure is that f3 is not flat. 



as shown in even if the intermittency range in the [0, 5] interval is very narrow. The 
second, is that the prediction (3 = D = .7565 is correct within two percent at a = and 
about five percent at a = 5. The third, and most important, that the scaling function j3 
is roughly invariant from case to case. 

These results are intriguing: the coincidence of the curves in Fig. 5 suggests that 
the spectrum of generalized dimensions Dq must play some role in determining /3(a): 
the fractal measures a - d seem to have little in common beyond having the same flat 
thermodynamics. Nevertheless, precisely because in these cases Dq is flat, neither eq. 
(P^D, nor any general relation of the kind /?(a) = Dq(^a)i with q an as yet unknown 
function of a can hold rigorously. 

Let us now come to a favourable example: we can construct a class of measures for 
which the renormalization eq. ( |3lD is exact these are the equilibrium measures of the 
Julia sets generated by the polynomials 

P{z) =z^ -\ (37) 

where A > 2 is a real constant. As we have already remarked, the Jacobi matrices 
for these problems can be constructed by a stable recursion algorithm [^. Non- 
linearity of the I.F.S. maps stemming from eq. (|37|) as inverse branches of P{z) can be 
treated by considering sufficiently high iterations P*-'\ and a theory perfectly analogous 
to (|30| - |35D can be carried out, with the same result. 

In fig. 6 we make the usual comparison between the moment scaling function, /?(a), 
and thermodynamics [^|, Di^a'- the curves coincide within numerical precision! There- 



fore, one can conclude that discrepancies from eq. (|35|) are due to the non-exactness of 
eq. ( pTD when a spectral measure is approximated by I.F.S., except for the case of Julia 
measures, which are known to have strong algebraic properties. 

Incidentally, we note that for Julia sets, the invariant measure coincides with the 
measure of the asymptotic distribution of the zeros of the associated orthogonal polyno- 
mials, the latter being also the physicists' global density of states. Might it be that the 
correct quantity entering eq. (|35D is this second measure ? The analysis of the I.F.S. 
data presented here seems to exclude this case, although we cannot exclude that this role 
is played by yet another spectral measure still to be determined. 

7 Conclusions 

We have presented a stable algorithm for the determination of lattice Hamiltonian op- 
erators possessing a given spectral measure, in the class of linear I.F.S. This algorithm 
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consists of a recursive determination of the associated Jacobi matrix, in the framework 
of the theory of orthogonal polynomials. 

The Hamiltonian operators determined in this way are characterized by almost peri- 
odic coefficients: since I.F.S. measures approximate arbitrary well any measure supported 
on a Cantor set, this fact might lead to a proof that almost periodicity is always associated 
with this kind of spectra. 

In a quantum mechanical context, the Jacobi matrices studied here can be employed 
as models of almost-periodic systems: the dynamical properties of such systems can be 
studied in their essence, having extracted the crucial information on the related spectral 
measures. We have shown that connections between spectral properties and dynamics go 
far beyond the conventional RAGE theorems: in particular, delocalization of particle's 
position along the lattice basis can be described by a scaling function j3 governing the 
moments of order a of the position operator. Non-constancy of this function translates 
mathematically the phenomenon of quantum interference. 

We have derived an intriguing relation, (3{a) = -Di-a, linking dynamics and the 
thermodynamical properties of the spectral measure: considering the Jacobi matrices 
associated with Julia sets we have constructed a family of quantum systems for which 
the relation is exact, and we have discussed the reasons for the discrepancies present in 
the general case. We believe that a further refinement of the results presented in this 
paper will lead to a profound understanding of the mathematical and physical properties 
of almost-periodic quantum systems. 

Finally, we remark that the Jacobi Hamiltonians considered in this paper are not sim- 
ple exotic curiosities, but can also describe time-resolved energy absorption in externally 
perturbed quantum systems, as well as electron dynamics in solid-state eterostructures 
like super-lattices [0, where by varying an alloy concentration along a deposition axis 
different spectral structures can be found [^. Here, our results may become relevant in 
several problems, like -for instance- the design of lasers and radiation detectors. 
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Figure Captions 



Fig. 1. 

Orthogonal polynomial ps{x) of the I.F.S. measure with maps {Si,f3i,Tri) — (|,0,|), 
(|, |, |), with a finite-resolution representation of the support of the measure obtained 
by plotting a large number of points on the attractor. Because of the finite size of points, 
this latter appears as a sequence of dashes. Only the symmetrical half is shown. In the 
inset, the beginning of the sequence of r„. The vertical scale ranges from zero to |. Lines 
are merely to guide the eye. 
Fig. 2. 

Discrete Fourier transform of the r„ sequence (n = 1, . . . , 2^^) associated with the I.F.S. 
of Fig. 1. The constant and tt frequencies exceed the vertical scale, and are not reported. 
Fig. 3. 

Plot of the phase $ of the discrete Fourier transform of the sequence associated with 
the I.F.S. of Fig. 1 and 2, to show the tt discontinuity close to the value of the main peak 
of Fig. 2. 
Fig. 4. 

Multi-fractal dimensions Dq (full diamonds) and Di (full squares) and dynamical expo- 
nents /3(0) (open squares) and /3(1) (open diamonds) versus contraction rate 5i, for the 
family of I.F.S. described in the text. 
Fig. 5. 

Scaling functions (3{a) for the four I.F.S. 's a - d described in the text: a: circles; b: 
squares; c: triangles; d: diamonds. 
Fig. 6. 

Scaling function P{a) for the Julia set measure with A = 2.2 (diamonds) and thermody- 
namical dimensions i^i-a (crosses). 
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